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ABSTRACT 

We present an X-ray study of the three most metal-deficient blue compact dwarf (BCD) galaxies known in the 
local Universe, based on deep Chandra observations of SBS 0335-052 (0.025Z Q ), SBS 0335-052W (O.O2Z ) 
and I Zw 18 (0.02Z Q ). All three are detected, with more than 90% of their X-ray emission arising from point- 
like sources. The 0.5-10.0 keV luminosities of these point sources are in the range (1.3-8.5)x 10 39 erg s . We 
interpret them to be single or a collection of high-mass X-ray binaries, the luminosities of which may have been 
enhanced by the low metallicity of the gas. There are hints of faint extended diffuse X-ray emission in both 
SBS 0335-052 and I Zw 18, probably associated with the superbubbles visible in both BCDs. The spectrum 
of I Zw 18 shows a O VIII hydrogen-like emission line. The best spectral fit gives an O overabundance of the 
gas in the X-ray point source by a factor of ~ 7 with respect to the Sun, or a factor of ~ 350 with respect to the 
O abundance determined for the HII region. 

Subject headings: galaxies: individual: IZw 18 — galaxies: individual: SBS 0335-052 — galaxies: individ- 
ual: SBS 0335-052W — galaxies: starburst — X-rays: galaxies 



1. INTRODUCTION 

Galaxy formation is one of the most fundamental prob- 
lems in astrophysics. To understand how galaxies form, we 
need to unravel how stars form from the primordial gas and 
how the first stars interact with their surrounding environ- 
ments. As there are no heavy elements in the early uni- 
verse, the thermodynamic behavior of the gas is essentially 
controlled by H2 cooling, and the first P opulatio n III stars 
are expected to be very massive (e.g.. | A.bel. Brvan. & Normanl 
l2002HBromm. Coppi. & Larsorl2002l) . While much progress 
has been made in finding large populations of galaxies at high 
(z > 3) redshifts (e.g., Steidel et al. 1996), truly young galax- 
ies in the process of forming remain elusive in the distant uni- 
verse. The spectra of those far-away galaxies generally indi- 
cate the presence of a substantial amount of heavy elements, 
implying previous star formation and metal enrichment. In- 
stead of focussing on high-redshift galaxies, another approach 
is to study massive star formation and its interaction with the 
ambient interstellar medium (ISM) in a class of nearby metal- 
deficient dwarf galaxies called Blue Compact Dwarf (BCD) 
galaxies, some of which are thought to be undergoing their 
first episode of star formation. 

The formation of the hot gas phase is one of the most funda- 
mental processes operating in the early period of galaxy for- 
mation and it has important consequences on the subsequent 
dynamical state and photometric evolution of the dwarf sys- 
tem. The generation of large amounts of hot (a few x 10 6 K), 
rarefied gas is the result of injection of energy and momentum 
into the cold ambient ISM by stellar winds from massive stars 
and supernovae. The starburst activity which fuels the ISM of 
young dwarf galaxies with hot X-ray emitting gas lasts about 
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10 million years. It is then followed by a long (> 1-2 Gyr) 
quiescent period of passive photometric evolution. Moreover, 
expansion of the hot ISM on scales comparable to the galac- 
tic scale length can lead to a rapid funneling of gaseous mass 
into the cold gaseous halo. The details of such a process are 
expected to be sensitively related to the energy injection rate 
into the ISM and the g eometry and robustness of the ambient 
cold gaseous medium (De Young & Heckmanl l 19941) . Exten- 
sive mass loss can lead to an expansion of the stellar popu- 
lation and to a morphological evolution of the dwarf galaxy 
(for example a dwarf irregular galaxy may lose its gas and 
evolv e into a dwarf elliptical) or even di srupt the system en- 
tirely jHillsll980):rYoshii & Arimotoll987l) . 

Here we present a study at X-ray energies of the three most 
metal-deficient star-forming dwarf galaxie s known in the 
local u niverse, the BCDs S BS 0335-052 (iThuan & Izotov 
1997 1 : llzotov et all Il997t IThuan. Izotov. & Lipovetskv 



1997 hereafter TIL), SBS 0335-052W ( Lipovetskv et al. 
1999T). and I Zw 18 ( see iPapaderos et alJ 120021 
iHunt. Thuan. & Izotovl 12003a! and references therein). 
X-ray emission provides an efficient and direct probe of the 
hot gas component, and can provide valuable constraints on 
the nature of the X-ray binary population. The number of 
high mass X-ray binaries (HMXBs) put limits on the BCD's 
supernova rate while the number of low mass X-ray binaries 
(LMXBs) contrains the total mass of forming stars. The 
relative proximity of these three BCDs [D = 54.3 Mpc for 
SBS 0335- 052 and SBS 0335-052W (TIL); D = 12.6 Mpc 
for I Zw 18 (Ostlin 2000)] , their extre mely low m etallicities 
[O.O25Z for SBS 03 35-052 Czotov et alJll999l) . 0.020Z Q 
for SBS 0335- 052W (iLipovetskv et alJU999h . and O.O2OZ 
for I Zw 18 (Izotov et al. 1999)], and high neutral gas 
content (see Pustil nik et alJ 120011 for SBS 0335-052 and 
SBS 0335-052W and Ivan Zeeetai] Il998l for I Zw 18) 
make them the best local approximations to primordial 
young galaxies. Their relatively bright apparent magnitudes 
(m B = 17.0 mag for SBS 0335-052, m B = 19.4 mag for 
SBS 0335-052W and m B = 15.8 mag for I Zw 18) make 
them much easier to study than the very faint and small 
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building-block dwarf galaxies at high-redshift. HI VLA 
mapping (Pusti lnik et all 120011) shows SBS 0335-052 and 
SBS 0335-052W to be embedded in a common very large 
HI envelope with an overall size of ~ 66x22 kpc. There 
are two prominent HI peaks separated in the east-west 
direction by 22 kpc (84")- The eastern peak is associated 
with SBS 0335-052, while the western peak is about a factor 
of 1.3 brighter in the HI line and is associated with the fainter 
S BS 0335-052W Like wise, HI VLA mapping of I Zw 18 
bv Ivan Zee et alJ 11998) shows that it is also embedded in an 
extended HI envelope. 

Optical HST imaging shows that nearly all star formation 
in SBS 0335-052 occurs in six very blue super-star clusters 
(SSCs) with absolute V luminosities between -1 1.7 and -14.7 
mag, and surface brightnesses > 100 times those of clusters 
and associations in normal HII regions. The SSCs are con- 
fined to a region of about 520 pc in diameter (TIL). There is 
a systematic color gradient from the brightest and bluest SSC 
at one end to the faintest and reddest SSC at the other end. 
Such behavior can be attributed partly to variable dust extinc- 
tion. Dust is clearly present, even in such a metal-poor en- 
vironment, as evide nced by the strong mid-infrared e mission 
of SBS 0335-052 JThuan. Savage. & Maddenl[T99l . How- 
ever, most of the color gradient is due to a systematic aging 
of the SSCs resulting from sequential propagating star forma- 
tion (with a velocity < 20 km s" 1 ) fro m the re ddest and oldest 
SSC to the bluest and youngest one (Panaderos et al. 1998) . 
From t he V-I colors and the models of Leitherer & Heckman 
( 1995), the age of the youngest SSC is ~ 4 Myr while that 
of the oldest SSC is ~ 30 Myr. On larger scales, the light 
of SBS 0335-052 is dominated by a patchy and filamentary 
very blue low surface brightness (LSB) emission. Spectro- 
scopic studies (Izotov et al. 1997) have revealed that ionized 
gas contributes ~ 30% of the emission from the LSB com- 
ponent, out to scales of ~ 3 kpc. The remaining ~ 70% of 
the light comes from a n underlying stellar population which 
Panad eros et alJ J1998I) have shown is not older than ~ 500 
Myr. A striking signature of the collective action of massive 
stars and SNe on the ISM of SBS 0335-052 is a large super- 
shell with radius ~ 380 pc to the North of the SSCs. TIL 
estimate the star formation rate in SBS 0335-052 from its Ha 
luminosity to be ~ 0.4 M Q yr" 1 . 

SBS 0335-0 5 2W h as been studied in detail by 
iLipovetskv et alJ d!999h . It has the same redshift as 
SBS 0335-052 and consists of at least three stellar clusters. It 
has also a very low metallicity of ~ 2% solar and is about 2.4 
magnitudes fainter than SBS 0335-052. Its very blue R-I 
colors also suggest a very young age of less than ~ 500 Myr. 

I Zw 18 is composed of two bright knots of star forma- 
tion separated by 5" and referred to as the brighter northwest 
(NW) and fainter southeast (SE) components ( Pan aderos et al.l 
120021) . Star formation proceeds differently in I Zw 18 as com- 
pared to SBS 0335-052. There are no SSCs in I Zw 18 and 
the derived SFR is one order of magni tude smaller (~ 0.04 
M^ yr" 1 ) than in SBS 0335-05 2.|papaderos eTall J2002I) and 
iHunt. Thuan. & Izotovl l l2003al) . using deep optical and near- 
infrared imaging, have put an upper limit of ~ 500 Myr for 
the age of I Zw 18. 

In f|2]we describe the properties of the X-ray sources and 
their optical properties. In f|3]we summarize our findings and 
discuss the X-ray properties of SBS 0335-052 and I Zw 18 
in the broader context of star formation in a low-metallicity 
environment. 



2. OBSERVATIONS AND DATA REDUCTION 
2.1. X-ray Observations 

SBS 0335-052 was observed on 2000 September 7 with 
ACIS-I (the imagin g array of the Advanced CCD Imaging 
Spectrometer; Garm ire et alJ20 02) aboard the Chandra satel- 
lite. The field of view of 16. '9 x 16. '9 was large enough to 
include also SBS 0335-052W. Events were telemetered in 
Very Faint mode, and the CCD temperature was -120° C. To 
compare the X-ray properties of SBS 0335-052 with those 
of another extremely metal-deficient BCD, we have retrieved 
from the Chandra archive the data for the BCD I Zw 18. This 
galaxy was observed on 2000 February 8 (P.I: D.J. Bomans) 
with ACIS-S (chips 235678 on, 8'5 x 8. '5 per chip). Events 
were telemetered in Faint mode, and the CCD temperature 
was -120° C. The ACIS pixel size is 0"492 while the half- 
energy radius of an on-axis point source is w 0"42. 

Analysis was done using primarily the CIAO V2.3 soft- 
ware provided by the Chandra X-ray Center (CXC), but 
also with FTOOLS and custom software. We removed the 
0."5 pixel randomization, cleaned the ACIS background us- 
ing Very Faint mode screening for SBS 0335-052, performed 
standard ASCA grade selection, and excluded bad pixels and 
columns. The total net exposure times were 61,042 s for 
SBS 0335-052 and 40,776 s for I Zw 18. The average qui- 
escent background rates for SBS 0335-052 and I Zw 18 were 
found to be 0.3 and 0.96 counts s" 1 per chip, respectively, 
in good agreement with ACIS-I and ACIS-S calibration mea- 
surements. The observation of SBS 0335-052 contained four 
background "flares", the intensity of each being about three 
times the quiescent level and lasting for 5-10% of the total 
exposure. The observation of I Zw 18 contained one flare 
about ten times the quiescent level lasting for 20% of the total 
exposure 5 . None of the flares were strong enough to war- 
rant removal. The resulting average background level is thus 
k 0.02 counts pixel" 1 and « 0.06 counts pixel" 1 , respectively. 

2.2. Astrometry of optical and X-ray images 

To compare the X-ray emission with the optical emis- 
sion, we have overlayed the X-ray contours over the V- 
band (F569W) HST WFPC2 image of SBS 0335-052 ob- 
tained by TIL and the V-band (F555W) HST WFPC2 image 
of I Zw 18 retr i eved f rom the archives and reprocessed by 
Panad eros etall J2002I) . As SBS 0335-052W does not he 
within the field of view of any HST image, we have used 
for optical comparison the CCD fi-band image obtained by 
Papaderos et al. ( 1998) with the 2.2 m Calar Alto telescope. 
The Calar Alto image has a large field of view (~ 17' x 17'), 
and contains enough stars for aligning the HST and X-ray 
images to the same astrometric reference frame. A similar 
alignment was performed for I Zw 1 8 us ing the CCD Z?-band 
image obtained by Papaderos et al. (2002), also with the 2.2 
m Calar Alto telescope. To align the optical images, sources 
were first extracted from the Calar Alto images using Sextrac- 
tor (v2.2.1: lBertin & Arnout3l 996). Optical sources from the 
Calar Alto images were matched to the GSC2.2 catalog for 
absolute astrometric alignment, giving a coincidence of 28 
sources within a 3(7 radius of 0."57 and a rms scatter of 0."25 
for SBS 0335-052, and 25 sources within a 3er radius of 0."62 
and a rms scatter of 0."23 for I Zw 18. In the same manner, 
optical sources from t he Calar Alto image w ere used to align 
the HST images (TIL; iPaoaderos etafll2 002 ) to the same ref- 
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erence frame; 51 sources were matched within a 0"82 radius 
and arms scatter of 0."31 for SBS 0335-052, while 21 sources 
were matched within a 0."49 radius and a rms scatter of 0"19 
forlZw 18. 

To align the Chandra image, we first extracted the X-ray 
sources using the CI AO wavelet algorithm WAVDETECT on 
the 0.5-8.0 keV imag es with a probability threshold of 10~ 7 
dFreeman et alJ l2002). which is equivalent to the detection of 
approximately one false source over the entire ACIS array. 
Eighty sources were detected in total in the SBS 0335-052 
field. Positions and counts for the sources were determined 
using fixed circular apertures with 95% encircled energy radii 
(based on point spread functions (PSF) derived from MKPSF 
and the PS F library). Sixty - four s ources lie within the region 
imaged by Papaderos et all i 19981) around SBS 0335-052, of 
which 18 X-ray sources have plausible optical counterparts 
within a 3a radius of 0."9, including both SBS 0335-052 and 
SBS 0335-052W. Less than 2 false matches are expected. 
Aligning the X-ray image to the optical astrometric reference 
frame based on these SBS 0335-052 X-ray/optical matches, 
we find that the astrometric accuracy of the X-ray positions 
has a rms scatter of 0"42. Likewise, 16 s ources lie within the 
region imaged by Papaderos et al. (2002) around I Zw 18, but 
only four X-ray sources have plausible optical counterparts 
within a 3a radius of l."7, including I Zw 18. This paucity 
of matches is likely due to the relatively shallow depth of the 
B-band Calar Alto image. Less than 1 false match is expected. 
Aligning the X-ray image to the optical astrometric reference 
frame based on these I Zw 18 X-ray /optical matches, we find 
that the astrometric accuracy of the X-ray positions has a rms 
scatter of 0."55. 

2.3. X-ray Analysis 

Based on our astrometric solutions, we present X-ray con- 
tours of SBS 0335-052, SBS 0335-052W, and I Zw 18 in 
Figure^ overlaid on optical images. SBS 0335-052 consists 
of a faint X-ray point source (29 counts) lying «0."3 north of 
SSC 2. We follow here the SSC numbering scheme of TIL. 
Because the X-ray position is no more accurate than 0."42, 
the data is consistent with the X-ray source being physically 
associated with SSC 2. The elongated shapes of the faint X- 
ray contours to the north suggest that there may also be X-ray 
emission associated with SSC 3 and the pair of SSC 4 and 5. 
There is also a very faint X-ray source to the northeast which 
appears to be associated with the supernova cavity seen in the 
optical image (TIL). Excluding the source near SSC 2, the 



total X-ray emission of the fainter features is <23 counts (Ta- 
ble 1), making it virtually impossible to distinguish between 
point-like and truly extended emission. 

Although SBS 0335-052W is 2.4 mag fainter in B than 
SBS 0335-052, it is 2.8 times brighter in the X-ray. It con- 
sists of two faint X-ray point sources. The fainter one appears 
to be spatially coincident with the brightest of the two optical 
clusters in SBS 0335-052W. The brighter X-ray source lies 
« 1 ."8 to the northeast and does not appear to have an optical 
counterpart. 

I Zw 18 is only 1.2 mag brighter in B than SBS 0335-052, 
but its X-ray flux is more than 10 times higher than that of 
SBS 0335-052. It consists of a single X-ray point source 
(471 counts) and a very faint extended X-ray component (22 
counts). The point source is spatially coincident with the 
bright NW component of I Zw 18. The faint extended compo- 
nent appears to lie along the main body of I Zw 1 8 with some 
X-ray emission centered on the northwest component and 
some in between the n orthwest and southeast optical peaks. 
Bomans & Weis ( 2002) in a preliminary reduction of the same 
Chandra data also found the X-ray emission of I Zw 18 to be 
composed of a strong point source and fainter extended emis- 
sion, although we fail to recover the SW and NE extensions 
these authors claim to see. 

X-ray variability is an extremely useful method for con- 
straining the size of the emitting region. Because Chan- 
dra records the position, energy, and time of all incoming 
photons, we can use this information to search for short- 
term variability over the duration of our observation. We 
searched for variability from the different X-ray components 
of SBS 0335-052, SBS 0335-052W, and I Zw 18 using the 
Kolomogrov-Smirnov (KS) test. While no significant vari- 
ability was found, the limited statistics do not strongly con- 
strain the temporal nature of the sources. 

X-ray spectra can provide information about the mecha- 
nisms by which the X-ray emission is gener ated. The X-ra y 
spectra were analyzed using XSPEC (vl 1.2 lArnaud l 1996). 
Unless stated otherwise, spectral parameter errors are for the 
90% confidence level, assuming one parameter of interest. 
The X-ray fluxes and absorption-corrected luminosities for 
all sources were calculated from spectral fitting using XSPEC. 
Importantly, none of the sources of interest suffer from pileup. 

To constrain the fit of the X-ray spectra, it is important 
to use accurate HI column densities for each of the objects. 
The HI c olumn density toward s SBS 0335-052 has been de- 
rived by ( Thua n"& Izotovin"997h through fitting the wings of 
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Table 1. X-ray Emission from SBS 0335-052, SBS 0335-052W, and I Zw 18 



(1) 




(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 


(10) 


Source 




Position 


Counts 


Model 


N» 


r/kT 


Fit/DOF 


Fx 


Lx 


Comments 


SBS 0335- 


-052 


033744.1-050239.5 


29.3 ± 6.5 


POW 


6.8 (< 16.3) 


2 1 +L5 


24.8/24 


6.1 


3.5 


Point Src. 










RAY 


5 9+ 63 


3.6 (> 1.2) 


24.7/24 


5.2 


2.8 












POW 


7.0 (fixed) 


2 2+06 


24.8/25 


5.7 


3.5 












RAY 


7.0 (fixed) 


L - '-1.3 


24.6/25 


4.5 


2.8 








033744. 1-050239.5B 


8.4 ±5.0 


RAY 


7.0 (fixed) 


1.0 (fixed) 




0.6 


0.64 


Extended 


SBS 0335- 


-052W 


033738.5-050236.5 


82.4 ± 10.2 


POW 


5 2 + ' i J 
J.A-2.7 


z - S -0.8 


41.1/56 


10.3 


8.5 


Point Src. #1 










RAY 


3 1+ 2 - 3 


2 o+ 2 - 2 


41.6/56 


9.6 


5.2 








033738.4-050237.3 


36.4 ±7.1 


POW 


2.3 « 7.1) 


1 9+' 1 


21.9/30 


6.3 


2.8 


Point Src. #2 










RAY 


1.3 (< 3.0) 


5.4 (> 1.9) 


22.0/30 


5.9 


2.4 




IZw 18 




093401.9+551428.4A 
093401. 9+55 1428.4B 


469.5 ± 21.7 
22.9 ± 6.9 


POW 
RAY 
VRAY 
RAY 


i 44+0.38 
0'.87|1 
0.94«| 
1.31 (fixed) 


2.01+OM 

4 06+ l|) ™ 
—1.19 

4 28+ 2 - 25 
1.0 (fixed) 


18.1/20* 
23.0/20* 
8.1/19* 


72.1 
66.6 
70.4 
2.0 


1.6 
1.4 

1.5 
0.053 


Point Src, 0.65 keV line? 

Z° = 7.0 + i 2 3 2 Z° 
Extended 



NOTE. — Column 1: Source name. Column 2: Source position given as CXOU JHHMMSS.S+DDMMSS.S. Column 3: Background-subtracted 0.5-8.0 keV counts accumulated 
over 60. 1 ks (SBS 0335—052) and 40.8 ks (I Zw 18). Aperture photometry was performed using 95% encircled-energy radii for 1.5 keV for point sources, and individual background 
regions were selected adjacent to each source as noted in ^2] The standard d eviation for th e source and background counts are computed following the method of Gehrels ' 1986) and 
are then combined following the "numerical method" described in § 1.7.3 of Lvons ' 1991 ). Column 4: Spectral model used to fit data. POW indicates an absorbed power-law model, 
whereas RAY (VRAY) indicates an absorbed Raymond-Smith thermal plasma model (with variable O abundance) ' Raymond & Smith 1977). Columns 5 and 6: Neutral hydrogen 
absorption column density (M0 in units of 10 21 cirT 2 . Photon index (r) or thermal plasma temperature {kT in units of keV) as determined from the best-fit absorbed power-law or 
thermal plasma models to the AOS spectra. Also listed are the 90% confidence errors calculated for one parameter of interest (Ax 2 = 2.7). Column 7: Goodness of fit/degree of 
freedom. For SBS 0335-052, fitting was performed with the c-statistic, while for I Zw 18 the \ 2 statistic was used (denoted by "*"). Columns 8 and 9: Observed 0.5-10.0 keV fluxes 
in units of 10~ 15 erg em" 2 s~ ! and absorption-corrected 0.5-10.0 keV luminosities in units of 10 39 erg s™ , assuming the best-fit model parameters given in Columns 5 and 6. Column 
10: Comments. 



the damped La profile in a HST/GHRS UV spectrum. They 
found N m = (7.0 ± 0.5) x 10 21 cm" 2 . This is nearly 10 times 
larger t han the peak HI col umn density of 7.4xl0 20 cm" 2 
which (Pustilni k et alJl2001l) found in their VIA HI map of 
SBS 0335-052. This large difference can be undertood as 
a beam-smearing effect in the VLA map. The synthesized 
VLA beam (20"5 x 15"0) is considerably larger than the 
GHRS aperture (2"0 x 2."0), and column densities in the di- 
rection of structures smaller than the VLA beam will be artifi- 
cially diminished. Thus, for the X-ray point source associated 
with SBS 0335-052, the HI column density obtained from the 
spectrum through the small GHRS aperture is the appropri- 
ate one to use. No comparable HST/GHRS spectrum exists 
for SBS 0335-052W. The VLA map (Pustilnik eTal]l200l 
shows that the HI column density towards SBS 0335-052W 
to be 1.35 larger as compared to SBS 0335-052. For lack of 
better information, we have simply scaled the GHRS HI col- 
umn density towards SBS 0335-052 by that factor to obtain 
the HI column density towards SBS 0335-052W. However, 
we should remember that, because of the large VLA beam, 
we do not know the exact value of N(HI) in the inner 2" of 
SBS 0335-052W, where the X-ray emission is located. For 
the HI column density towards I Zw 18, we have adopted the 
value N m = (3.5 ±0.5) x 10 21 cm" 2 derived by (Ku nthet all 
119941) in the same way as for SBS 0335-052, by fitting 
the wings of the damped La profile in a HST/GHRS spec- 
trum of the NW component, obtained through a 2."0 x 2"0 
aperture. This value is en tirely consistent with the value of 
2.1 x 10 21 cm" 2 derived by Vidal-Madjar (2000) from fitting 
the Ly/3 absorption profile, with that of (2.0±0.5) x 10 21 cm" 2 
obtained by Lecavelier des Etangs et al. (2003) and that of 
2.2+^ x 10 21 cm" 2 obtained by Aloisi et al. (2003) by fitting 
the blue wing of the Ly/3 line and the profiles of several other 
lines of the HI Lyman series. The last three determinations 
are base d on FUSE spectra o btained through a 30."0 x 30."0 
aperture. [Brown et al.l J2002I) have obtained STIS data of I Zw 



18 at the very high spatial resolution of 0."5. They discov- 
ered significant inhomogeneity in the HI gas, with a peak as 
high as Nm ~2 x 10 22 cm" 2 . However, this is not the appro- 
priate value to use here as this HI peak is associated with the 
SE component, while th e X-ray source is a ssociated with the 
NW component which JKunth et aUll994l) observed. VLA 
HI maps of I Zw 18 have been obtained by van Zee et al. 
(1998). They also found a peak HI column density at the SE 
component of ~ 3 x 10 21 cm" 2 . Again, because the highest 
resolution map of van Zee et al. (1998) has a beam size of 
5 "2 x 4."8, beam smearing effects have decreased the real HI 
peak column density by a factor of nearly 10. 

The ACIS spectra of the sources in SBS 0335-052, 
SBS 0335-052W, and I Zw 18 were fit with two types of 
models, absorb ed therma l plasma models {zwabs+raymond; 
Ravmond & Smith 1977) and absorbed power-law mod- 
els {zwabs+pow). Given the limited counting statistics, 
we performed the spectral analysis using the unbinned, 
background-subtracted source spectra and the Cash statistic 
jjC^gh|1979), w hich is well suited to low-count sources (e.g., 
Nousek & Shue 1989). Note that we expect < 2 background 
counts in the source aperture, so spurious residuals from spa- 
tially varying background should be negligible. One limita- 
tion of the Cash statistic is that it does not provide a reliable 
quality-of-fit criterion (like the \ 2 statistic) to compare differ- 
ent models. Therefore, we established the quality of the spec- 
tral results via visual inspection using the binned spectrum. 
None of the sources have enough counts to constrain strongly 
both the column density and the temperature or photon index. 
In our spectral fits using a thermal plasma model, we have 
assumed an abundance of Z = 0.025Z Q for SBS 0335-052, 
and Z = O.O2Z for SBS 0335-052W and I Zw 18, in ac- 
cord ance with the ab undances determined spectroscopically 
(see llzotov et al.|fl999l for SBS 0335-052 and I Zw 18, and 
iLipovetskv et all 19991 for SBS 0335-052W). 

The spectrum of the point source in SBS 0335-052 has only 
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«29 counts and is not strongly constrained by either the ab- 
sorbed power-law or thermal plasma models when the column 
density Nu is left as a free parameter in the fit (see Table 
If we fix A^h to the H I value of 7 x 10 21 cm 2 obtained from fit- 
ting the Lya profile (iThuan & Izotovlll997l) . we find that the 
spectrum is moderately soft, with best-fit values of the pho- 
ton index T = 2.2^jg (V = a + 1, where a is the slope of the 
power-law fit ) or kT = 2.7+[ 6 3 ' 6 keV. These values are con- 
sist ent with those of b lack-hol e X-ray binarie s (e.g., T ~ 1.7- 
2.2; Foschini et al. 2002; Roberts et al. 2003), but are not con- 
sistent with a hot gas component heated by supernovae which 
has typically kT < 1 keV. The point source has an unabsorbed 
X-ray luminosity of (2.8-3.5) x 10 39 erg s" 1 . The spectrum 
and lack of variability do not provide further constraints on 
the nature of this source. 

The spectrum of the brighter component point source (#1) 
in SBS 0335-052W is acceptably fit by both the absorbed 
power-law or thermal plasma models. The best fit param- 
eters are N m = (5.2! 3 ; 3 ) x 10 21 cm" 2 , T = 2.8!g| and N m = 
(3.1!?; 3 ,) x 10 21 cm" 2 , kT = 2.0^j keV, respectively (see Ta- 
ble Q). The A^h vs. T and A^h vs. kT confidence con- 
tours for SBS 0335-052W #1 are given in Figure indi- 
cating the relative accuracy of the derived parameters. The 
soft intrinsic spectrum is similar to that of the point source 
in SBS 0335-052, and is not consistent with the typical 
hot gas component seen in other galaxies. The best fit val- 
ues of A^h in both models are higher than the beam-diluted 
Nm = 1 .0 x 10 21 cm" 2 obtained from the VLA HI map. How- 
ever, they are lower than the value of 9.5 xlO 21 cm" 2 ob- 
tained from scaling the Nm obtained from the fit of the 
damped La profile in SBS 0335-052 by the ratio of radio 
HI column densitie s of 1.35 between SBS 0335-052 and 
SBS 0335-052W dPustilnik et al.ll200l . The point source in 
SBS 0335-052W has an unabsorbed X-ray luminosity in the 
range (5.2-8.5) xlO 39 ergs" 1 . 

The spectrum of the fainter component point source (#2) 
in SBS 0335-052W is also acceptably fit by both the ab- 
sorbed power-law or thermal plasma models. The best fit 
parameters are N H = 2.3(< 7.1) x 10 21 cm 2 , T = 1.9^ and 
N H = 1 -3« 3.0) x 10 21 cm 2 , kT = 5.4(> 1 .9) keV, respectively 
(see Table Although the column density is not well con- 
strained, the upper limits are in good agreement with those 
from point source #1. The intrinsic spectrum of this source is 
the hardest of the three discussed so far, but it is still relatively 
soft. This point source has an unabsorbed X-ray luminosity in 
the range (2.4-2.8) x 10 39 erg s" 1 . 

The spectrum of I Zw 18 has much better photon statistics 
than SBS 0335-052, so we grouped the X-ray spectra in 20- 
count bins and fit the data using x 2 - The absorbed power-law 
model provides a slightly better fit to the continuum than the 
absorbed thermal plasma model, although both are formally 
acceptable. The best fit parameters are A^h = (1. 44^37) x 
10 21 cm" 2 , T = 2.01^ and N B = (0.87!^) x 10 21 cm" 2 , 
kT =4.06^[ jg keV, respectively (see Table [Q. Interestingly, 
we find large residuals below « 1 keV, particularly around 
w 0.65 ke V, which is most likely cau sed by one or more emis- 
sion lines. Bomans & Weis (2002) also found this line, ac- 
counting for about 3% of the total flux, which they attributed 
to the OVIII hydrogen-like emission line. 

One way to account for this particular emission feature is an 
O enhancement in the X-ray gas. If we model the spectrum 
with an absorbed thermal plasma model, varying the O abun- 
dance but keeping all other elements fixed to Z = 0.02Z Q , we 



find an improvement over the simple thermal model fit at the 
99.99% confidence level (Figure [3}- The best fit parameters 
are 

N H = (0.94^11) x 10 21 cm" 2 , kT = 4.28! 2 ^ keV, Z° = 
7.0^ 2 3 2 Z°, (see Table[l). 

To check for the possibility that a-elements other than O 
may also be contributing to the shape of the X-ray spectrum, 
we have also investigated a model where the abundances of 
the a-elements Ne, Si, S and Ar vary in concert with the O 
abundance, but in such a way that the a-element-to-oxygen 
abundance ratios Ne/O, Si/O, S/O and Ar/O remain constant, 
equal to the values derived bv llzotov & Thuanl(ll99 9) for low- 
metallicity BCDs (see their Table 6). This model gives a 
worse fit to the X-ray spectrum. We conclude that the O abun- 
dance of the X-ray gas in the point source in I Zw 18 is likely 
to be enhanced by a factor of ~ 350 with respect to the O 
abundance determined for the HII region gas, while the other 
a-elements have abundances that are consistent with no en- 
richment with respect to the abundances in the HII region. 

The point source in I Zw 18 has an unabsorbed X-ray lu- 
minosity between 1 .3-1 .5 x 10 39 erg s" 1 , a factor of 2-4 times 
smaller than those found for SBS 0335-052. In addition to 
the point sources discussed previously, both SBS 0335-052 
and I Zw 18 also show faint extended sources. The flux of 
the extended source in SBS 0335-052 is only about 10% of 
the 0.5-10 keV flux of the point source, and that of the one 
in I Zw 18 is only about 2.7% of the flux of the point source 
(Table 1). The physical nature of these extended sources will 
be discussed more in detail in the next section. 

3. RESULTS AND DISCUSSION 

Our main findings are the following: 

1) X-ray emission is detected from both SBS 0335-052 
and SBS 0335-052W, and from I Zw 18. The 0.5-10.0 keV 
luminosities of the three BCDs are in the range 1.3-8.5 x 
10 39 erg s" 1 . They are dominated by point sources, although 
at the distances of these objects and with the limited statis- 
tics, it is difficult to rule out emission from a collection of 
fainter point or compact extended sources. If these sources 
are single objects, their luminosities would place them in the 
rang e of the so-called ultra luminous X-ray sources (ULXs; 
e.g., iMakishima et al.l2000l) 

The X-ray spectra of the point sources are well fitted 
by moderately soft power-laws, typical of X-ray emission 
from HMXBs or ULXs. The high X-ray luminosities of 
these sources could be the result of the lower metallicities 
of the BCDs studied here (if accretion is via stellar winds), 
or beaming (as is thought to occur in some ULXs; e.g., 
King et al. 2001). The mean luminosity of galactic H MXBs is 
~ 5 x 10 36 erg s" 1 (Ivan Paradiis & McClintockl 19951) . The X- 
ray sources in SBS 0335-052, SBS 0335-052W and I Zw 18 
are respectively about 600, 1300 and 300 times brighter. For 
comparison, the HMXBs in the Magellanic Clouds (the LMC 
and the SMC have respectively 1/3 and 1/8 of the Sun's metal- 
licity) are about 50 times b righter than their galactic counter- 
parts (Ivan Paradiis & McClintockll9 95). 

In the event of X-ray emission produced by accretion 
onto a compact object by a stellar wind, there are two 
main reasons which may explain the boosting of the X- 
ray lu minositi es of HMXBs in low-metallicity systems 
(van Paradiis & McClintock 1995). First, the X-ray heating 
of the gas as it falls toward a compact object depends strongly 
on the atomic number Z, since the photoelectric cross-section 
goes as Z 4 . Thus in low-metallicity systems, the X-ray heat- 
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FIG. 2.— (Left) N H vs. T and (Right) N H vs. &r confidence contours for SBS 0335-052 West #1. Contours delineate the 66%, 95%, and 99% confidence 
regions. 
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FIG. 3. — (Left) X-ray spectrum of I Zw 18 point source, modeled with an absorbed power-law of slope T = 2.01 and Nh = 1.44 X 10 21 cm~ 2 . Note the large 
residuals around 0.65 keV suggesting the presence of an OVIII hydrogen-like emission line. (Right) Same spectrum, modeled with an absorbed thermal plasma 
of kT = 4.28 keV, Nu = 9.4 X 10 20 cnr 2 , and an O abundance of Z° = 7.0ZS, which significantly improves the fit in the 0.65 keV region. 



ing is less, resulting in a larger accretion rate because it is 
less impeded by heating. Second, the accretion rate dM/dt in- 
creases strongly as the stellar wind velocity v of an OB star 
as the orbit of the compact object decreases (dM/dt is propor- 
tional to v" 4 ). It is known that the terminal wind velocity in 
low-m etallicity systems drops significantly. Thu an & Izotovl 
(1997) have measured a terminal velocity of ~ 500 km s 
for SBS 0335-052, as compared to 2000-4000 km s" 1 for 
the Galaxy and the LMC, and 1200-1500 km s _I for the 
SMC. Thus the high X-ray luminosities of the sources in 
the metal-deficient galaxies SBS 0335-052, SBS 0335-052W 
and I Zw 18 (i.e., Lx > 10 39 erg s" 1 ) do not necessarily require 
extreme black hole masses. 

The X-ray luminosities found here are higher than 
the point-source luminosities observed in less metal-poor 
dwarf starburst dwarf galaxies such as NGC 3077 (Z Q , 
£m— a.nkBV ~ 2-5 x 10 38 > 15% from point sources; 
Ott. Martin . & Walterll2003l) . IC 10 (Z /4, L a ws n key « 2 x 
10 38 , > 95% from a single point source; Bauer & Brandt 
120031 NGC 4449 ( Z^/4, Lm-snfr.v ~ 2 x 10 39 , ~ 60% 
from point sources; ISummers et alJl20 03) and NGC 1569 
(Z^/5, Lm— finkpv ~ 1 x 10 39 , ~ 30% from point s ources; 
Martin, Kobulni ckv. & Heckmanl l2002t iHeike et alJ 120031) . 
They are comparable to the point source luminosity in Holm- 
berg II (Z^/1 0, L 1—ino k ,.v » 1 x 10 40 , > 75% fr om a single 
point source; Mivaji, Lehm ann. & Hasingerl2001l) . 



Interestingly, iHunt et all J2003bl) have detected a non- 
thermal radio source associated with SSCs 1 and 2. Its flux 
density is ~ 0.6 mJy at 5 GHz, corresponding to a non- 
thermal radio luminosity of ~ 2 x 10 20 W Hz" 1 . Within the 
uncertainties in the X-ray and radio positions, the non-thermal 
radio source can be associated with the X-ray point source. 
However the nature of such a radio+X-ray emitting object is 
not clear. While its X-ray luminosity is similar to that ex- 
pected from a microquasar associated with a medium-mass 
black hole, its radio luminosity is about a mill ion times g reater 
than that observed in such objects {Mirabel & Rodriguez 
1999). It can also be that the X-ray and radio sources are 
physically distinct objects, although they are both associated 
with the massive starburst episode that gave birth to the SSCs. 

We note also that the time scale for HMXB formation of 
3-10 Myr after the onset of the starburst is consistent with the 
ages of 5 Myr or y ounger of SSCs 1 and 2 in SBS 0335-052 
( Vanzi et al. 2000) a nd of ~ 3 Myr of t he NW component in 
I Zw 18 JHunt. Thuan. & Izoto\l2003'3) . 

2) In addition to the point source emission in SBS 0335- 
052 and I Zw 18, there appears also to be faint extended dif- 
fuse X-ray emission from hot gas, hints of which can be seen 
in the smoothed X-ray contours of Figure ^ particularly in 
I Zw 18. This emission, if truly diffuse, can be understood 
in both BCDs as produced by the hot gas filling the cavity 
of a superbubble carved out in the ISM of the BCD by stel- 
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lar winds and supernovae. Indeed, HST optical imaging by 
TIL shows a large superbubble cavity of radius ~ 380 pc in 
SBS 0335-052, northeast of SSC 3 (see their Figure 1), where 
the diffuse X -ray emission appears to lie. Using equations 
2, 4 and 6 of Martin & Kennicut ( 1995), the X-ray luminos- 
ity Lx of the superbubble is proportional to /j 1 - 65 n °- 66 7 ,2 - 66 i 
where R is the radius of the superbubble, n the ambient ISM 
density and T the gas temperature. Adopting R = 380 pc, 
n = 12 cm -3 and T = 7.6 x 10 6 K, the predicted L x is in 
good agreement with the observed diffuse X-ray luminosity 
of 6.4 x 10 38 ergs s" 1 . As for I Zw 18, the diffuse X-ray 
emission lie between the NW and SE centers of star f orma- 
tion. HST optical imaging by Hunte r* & Thronsonl(ll995l) sug- 
gests the presence of a supernova cavity there of radius 120 
pc (see their figure 3). With R = 120 pc, n = 12 cm -3 and 
T = 6.0 x 10 6 K, L x is equal to 5.3 x 10 37 ergs s" 1 , again in 
good agreement with the observed diffuse X-ray luminosity 
inIZwl8. 

3) The spectrum of I Zw 18 shows a stron g emis sion line at 
~ 0.65 KeV, identified bv lBomans & Weisl ll2002l) to be a O 
VIII hydrogen-like line. The best spectral fit (Figure [3} gives 
an overabundance of oxygen by a factor of ~ 7 with respect 
to the Sun and of ~ 350 with respect to the other heavy ele- 
ments. It is interesting that this line is so strong be cause FUSE 
observations of the HI interstellar gas in I Zw 1 8 (lAloisi et alJ 
l2003ULecavelier des Etangs et allM B) shows the ISM to be 
as or more metal-deficient than the HII gas, i.e. to be less than 
2% of the solar metallicity. The presence of a sizable amount 



of oxygen in the X-ray spectrum of the point source is most 
likely due to oxygen-enrichment by either a companion star 
(in the case of an X-ray binary) or undiluted supernova ejecta 
(in the case of hot gas). 

Further progress in our understanding of the X-ray emis- 
sion in SBS 0335-052, SBS 0335-052W and I Zw 18 will 
require substantially better photon statistics and higher spa- 
tial resolution. A long observation with XMM-Newton, or 
Constellation-XIXEUS in the future, would be able to provide 
better spectral and temporal constraints that might confirm the 
nature of the brightest point sources as true single compact ob- 
jects or a collection of X-ray binaries and hot gas. However, 
such observations will not have the necessary spatial resolu- 
tion for separating the different emission components seen by 
Chandra. Spatially resolving these objects with good sensi- 
tivity is beyond the capability of Chandra and must await the 
launch of Gene ration-X or a similar high-resolution X-ray ob- 
servatory. 

Putting our results in a cosmological context, we expect that 
the X-ray emission from high-reshift primordial galaxies to 
come mainly from very high luminosity HMXBs. 
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